The spherical phylon group and invariants of 

v6"- the Laplace transform. 

as : 
as . 

A.L. Carey 1 , M.G. Eastwood 1 , P.E. Jupp 2 and M.K. Murray 1 

> ■ 

February 7, 2008 

vo : 

(N ■ 

1 Department of Pure Mathematics, University of Adelaide, 
>. I Adelaide SA 5005, Australia 

2 School of Mathematical and Computational Sciences, 
q ■ University of St Andrews, North Haugh, St Andrews KY16 9SS, U.K. 

■ 

Abstract 

(3jT)j We introduce the spherical phylon group, a subgroup of the group 

^jjy of all formal diffeomorphisms of M d that fix the origin. The invari- 

"■^J , ant theory of the spherical phylon group is used to understand the 

invariants of the Laplace transform. 

& ; 1 Introduction 

There are various contexts in which it is of interest to find approximations 
for large n to integrals of the form 



/ 



<dyLY>{-nf(x)}b(x)dx l . . . dx d , (1.1) 



where / and b are smooth functions on W d with / having a Taylor series at 
of the form 

~^fij%^ 7^ fijk% %~*% + . . . , (h2) 

and 6(0) ^ 0. Here the matrix (f^) is required to be positive-definite and we 
assume throughout the Einstein summation convention. 
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Approximations to integrals such as ( |1 . 1| ) are useful in statistical infer- 
ence, see, for example Barndorff-Nielsen & Wood (1995) and Shun Sz McCul- 
lagh (1995). Laplace's method (see e.g. §6.4 of Bender & Orszag (1978) or 
§8.3 of Bleistein & Handelsman (1986)) provides an asymptotic expansion in 
n of (11.1ft of the form 

n- d / 2 {A (/, 6) +Ai(/,6K 1/2 + ...}, (1.3) 

The coefficients Aj depend on the function / and the measure bdx 1 . . . dx d . 
The behaviour of ( |1.3j ) under reparameterisation has been discussed in Pace 
& Salvan (1994) and Barndorff-Nielsen & Wood (1995). Since the function / 
and the measure bdx 1 . . . dx d do not depend on any choice of coordinates on 
M. d , the integral ( |1 . 1|) does not depend on any coordinate system. It follows 
that the coefficients Aj(/, b) in ( |1.3| ) are independent of the choice of coor- 
dinates on W 1 . In other words, the Aj(/, b) are diffeomorphism invariants of 
the pair (/, b). We show below that (if d > 1) these invariants are essentially 
the only invariants of the pair (f,b). By this we mean that if we have two 
such pairs (/, b) and (g, c) such that Aj(/, b) = Ai(g, c) for i = 0, 1, . . . then 
there is a diffeomorphism mapping (/, b) to (g,c). 

There are a number of settings in which we could attempt to prove this 
result. For example, we could take / and b to be defined globally on a d- 
dimensional manifold or locally in a neighbourhood of the origin in R d . We 
consider here the case that / and b are jets, or, because M. d has a canonical co- 
ordinate system, formal power series. Since the asymptotic expansion ( |1.3|) 
is obtained by Taylor expansion of / and b about 0, the coefficients Aj(/, b) 
depend on / and b only through their Taylor expansions at 0. Thus it is 
possible to define the Aj(/, b) when / and b are formal power series. This 
approach is consistent with Pace & Salvan (1994) and fits in with previous 
work, on the phylon group, of three of the authors contained in Barndorff- 
Nielsen et al. (1992), Carey & Murray (1990) and Carey et al. (1996). 
Expressions for the first few of the Aj are given in Pace & Salvan (1994) and, 
as might be expected, these rapidly become complicated as % increases. In 
the case that b = 1, a combinatorial formula for the Aj(/, b) can be obtained 
from equation (4) of Shun & McCullagh (1995). We also present a general 
formula for the Aj(/, b) in ( |2.1| ). We find that the Aj(/, b) are polynomials in 
the coefficients of / and b and rational in the square root of the determinant 
of the Hessian matrix of /. This is entirely analogous to classical invariant 
theory. 

We begin with the calculation of the general formula for the invariants 
and follow that with a proof of the theorem. 
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2 A general formula for the invariants 



First let us fix some notation. Given (/, b) let f(x) = f(x) — (l/2)f i jX l x : > and 
denote by the inverse of the Hessian matrix (fy). Then we can write 
the integral ( |1.1|) as 

f exp{—n-fijX l x : '}exp{—nf( y x)}b(x)dx 1 . . . dx d . 

Standard techniques from distribution theory (see Appendix ||) can be used 
to show that if % is odd then Aj(/, b) = and if % is even then 



(o~\d/2 C— iV 

M/ ' J ~ det(^)V2 2^ n r \ s ,V h^n-jJ , l^, 



where (l/ , r!)(/ / )i 1 ...j r is the coefficient of x 11 . . .x lr in the Taylor expansion 
of (/)', the /-th product of /, and k±, k^l ■ ■ ■ |Av+ s _i, k r+s runs through all 
partitions of i\, . . . , i s , ji, . . . , j s into subsets of size 2. For i = 0, 2 the formula 
(|2.1) for the invariants agrees with the results in Pace & Salvan (1994). 



3 The phylon group 

Denote by Ai(d) the set of all formal power series /, in d variables, of the 
form ( |1 . 2| ) with the matrix (fy) positive-definite and by V(d) the set of all 
formal power series b with 6(0) ^ 0. Denote by Ai k (d) the set of all formal 
power series / which begin in the same way as ( |1.2| ) but which are finite of 
degree k. There is a projection map 

/: M{d) -> M k (d), 

the 'fc-jet' map which just truncates a formal power series. Similarly we 
define V k (d) and j k : V(d) -> V k (d). Notice that j°(b) = b(0). 

At various places we shall need one non-trivial fact about formal power 
series. If h is a smooth function on a neighbourhood of in M. d then its Taylor 
series at defines a formal power series, which we denote by j°°(h). This 
map is well-known not to be injective. However it follows from a theorem of 
Borel (see e.g. p. 390 of Treves (1967)), that it is surjective, that is every 
formal power series is the Taylor series at of some smooth function. We 
shall use this fact several times. 

Denote by V(d) the set of all formal power series with values in M. d which 
have no constant term and an invertible linear term. An element ip in V(d) 
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has the form ip = . . . , ip d ) where 
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^( x ) = p jX 3 + j-^) k x j x k + ... 

with (ijjj) an invertible matrix. Under composition of formal power series, 
V(d) forms a group, called the phylon group. To prove this, one can either 
construct inverses explicitly using the fact that (ip*) is an invertible matrix 
or replace ip by a smooth function using Borel's theorem and apply the 
inverse function theorem. (For further details of V(d) and its representations 
see Barndorff-Nielsen et al. (1992), Carey & Murray (1990), Carey et al. 
(1996).) 

We can also define groups V k (d) of finite formal power series of degree k 
with no constant term and with invertible linear term. Their multiplication 
consists of composition followed by truncation to order k. With this multi- 
plication the jet map j k : V(d) — > V k (d) is a homomorphism. We denote the 
kernel of j k : V(d) — > V k (d) by V^(d). This is the normal subgroup of all ip 
with components of the form 

nx) = 5)x> + J^^.., k+1 * n ■ --^ k+1 + ■■■■ (3-1) 

Of particular interest is the first of these groups, V^ l \d). The quotient 
V(d)/V iyl \d) is just GL(d), the general linear group of all invertible d by d 
matrices. Because a linear transformation is a global diffeomorphism of M. d , 
the general linear group is also a subgroup of V(d) and it follows that V(d) 
is a semi-direct product of GL(d) and V^(d). 

The phylon group acts on the space A4(d) by ip sending / to ipf, where 

{^f){x) = f{r\x)) ) (3.2) 
and on V(d) by ip sending b to ipb, where 



(#)(x) = 6(^" 1 (x)) 



det( 



d(^- l )\ 
dx'j ' 



(3.3) 



Notice that a slight subtlety arises here in the interpretation of the modulus 
signs. If 

a = a(0) + ajX* + . . . 
is a power series with a(0) ^ then \a\ denotes the power series defined by 

, a if a(0) > 
-a it a(0) < 0. 
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The group V k (d) acts similarly on Ai k (d) and V k l {d). These group 
actions commute with the projections j k in the sense that 

3 k ^f)=3 k m\f) and = j k+1 &)j k (b). (3.5) 

We have noted in the introduction that if we define the invariants A; by 



the asymptotic expansion of ( |1 . 1|) then they are manifestly invariant under 
local diffeomorphisms of M. d fixing the origin. We need to know that when 
we consider the Aj as functions on J\4(d) x V(d) they are invariant under 
the action of V{d). This can be proved using Borel's theorem. We just 
replace the formal /, b and ip by functions /, b and a local diffeomorphism 
ip. Invariance follows from the fact that j°°{ipf) = if>f and j°°{ipb) = ipb. 
We wish to prove the following theorem. 

Theorem 3.1. Let d > 1. If (/, 6) and (g,c) are elements of Ai(d) x V(d) 
then 

Hf,b)=A t (g,c) for alii = 0,1,2,... 
if and only if there is a in V(d) such that 

(/,&) = (ipg,tj)c). 



We begin our proof of the theorem by noting that if / is a function, 
not just a formal power series, with Taylor series as in ( |1.2|) then Morse's 
Lemma (Milnor (1963)) shows that we can change co-ordinates so that / is 
the standard quadratic form 

q(x) = \\x\\ 2 . (3.6) 

By using Borel's theorem we can show that this result is true also for formal 
power series. That means that if / G J\4(d) then we can find a ip in V(d) 
such that ij)f = q. 

To understand when (/, b) and (g, c) can be transformed one into the 
other by the phylon group, we can begin by transforming them so that / and 



g both become q, where q is defined in equation (|3.6| ). Then we can try to 
transform b into c without changing q. That is, we act by the subgroup of 
the phylon group that fixes q. We will call this the spherical phylon group 
and denote it by S(d). The spherical phylon group is therefore the subgroup 
of all ip G V(d) such that 

*)|| 2 = IM| 2 . (3.7) 
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We can now see why the theorem is not going to work when d — 1. In 
that case the spherical phylon group is just 0(1) = {+1,-1}. We discuss 
the one-dimensional case in Section 

Before discussing the proof of Theorem |3J] we need to analyse the spher- 
ical phylon group in more detail. 

Notice that the intersection of GL(d) with S(d) is precisely the group 
0(d) of orthogonal transformations. Moreover, if ip, defined by 

il>(x) = i>)x j + -^) k x j x k + ..., 

is in the phylon group then it is easy to check that X = is orthogonal. 
Hence X -1 ^ is also in the spherical phylon group, so that S(d) is the semi- 
direct product of 0(d) and 5 (1) (d), where S^(d) = S(d) n V^(d). If we 
apply the condition for an element ip with components 

ij\x) = 6}x j + ^4>} k x j x k + ... 

to be in S(d), we see that this is so if 

fc+i 

2%v}.,.,, i2) + 5 i5 E4,.,v4 + ,., fc+2) = (3-8) 

r=2 

for all k = 1,2,..., where we use the convention that indices appearing 
between round brackets are symmetrised over. Define S^(d) = S(d) n 
V^(d). If ip is an element of S^ k '(d) then it satisfies ( |3.1| ) and so the condition 
in equation (|3.8|) reduces to 

S ^L. M = °- ( 3 - 9 ) 
Let us write Aj(6) for Aj(g, b). Again, straightforward calculations using 
moments of the multivariate normal distributions or direct application of 
( |2.1| ) shows that 

Hh) = * d/2 (-^I ct ^( & ) ; (3-10) 

where the i-th complete trace, ctr«, of b is defined by 

f if % is odd 

c r HJ - I J2 ki _ k , b klkl ... k]kj if i = 2] is even 1 ' 

and (l/r\)bi 1 ...i r is the coefficient of x 11 . . .x %r in the Taylor expansion of b. 
Theorem |3.1| is a consequence of the following theorem: 

Theorem 3.2. If b, c are elements of V(d) with d > 1 then there is an 
element of the spherical phylon group ip such that ip(b) = c if and only if 
Ai(b) = Ai(c) for all i. 
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4 Proof of theorem |St2 



In one direction the proof follows from the invariance of the Aj under the 
action of the phylon group. To prove the other direction, start with b and c 
in V(d) such that Aj(c) = Aj(6) for all i = 0, 1, . . . . We shall prove below 
the following Proposition: 

Proposition 4.1. If b and c are elements of V(d) such that Aj(c) = Aj(6) 
for all i and j k ~ 1 (c) = j k ~ l {b) then there exists ipk G S^ k '{d) such that 
J k (c) =j k (ip k b). 



Before proving Proposition |4.1| , let us see how it proves Theorem p.2| . It 
follows from ( |3.10| ) that c(0) = 6(0) so, by definition, j°(c) = j°(b). We now 
use Proposition |4.1| to manufacture a sequence of elements ipk G S^ k ' (d) for 
each k = 1, 2, . . . such that 

3 k (c) = j k (^ k . 1 ...^J 1 b). 

Note that when we multiply ipk-i ■ ■ ■ i>i by ipk on the left we have 

/(V^-i • • • = j k (^k)j k (ip k -i ■ ■ ■ V>0 
= /(V>*-i...Vi), 



because ipk G S^ k \d) = ker(j fc ). Hence multiplying ipk-i ■ ■ ■ ifti by iftk changes 
no terms of degree lower than k + 1. This means that there is a well-defined 
formal power series ip G V(d) such that for all k = 1, 2, . . . 

J fc W=J fc (^-i...^i). 

Then we have 

j fc (c) =j k {ipki) k -i---i>ib) 

= j k+1 (iPkTpk-i...<Pi)j k (b) 

=3 k+l m k {b) 
=/(#) 



for all = 1,2,.... Hence c = proving Theorem |3]2 
It remains to prove Proposition [4.1| . Write 

b=a+ B+b 
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and 

c = a + C + c, 

where a is a polynomial of degree of k — 1, B and C are homogeneous of 
degree k and j k {b) = = j fe (c). 
For any ■?/> in S^ k \d), write 

^(x) = ^ + j^y^.., k+ ^ 31 ■ ■ ■ * Jk+1 + o(M fe+2 ). 



Then 



so that 



dib l 1 

-^- = 5) + , . . . z 3 '* + 0(|x| fe+1 ) 



det (|^) = 1 + yr..,r' J ■ ■ ■ ^ + 0(1*0- (4.1) 



Note that because the power series in Qll|) is 1 at 0, its absolute value is 
itself (by ( |3.4[ )). Then from ( |3.3| ) we have 

dtp 1 



ip 1 b = {b o if)) det 
so that 



dx j 



= a (x) + B(x) + h -^-^\ h , ... jk x jl ■ ■ ■ % jk + 0(\x 



If we can find a 1 satisfying (|3.9| ) such that 

C(x) = B(x) + 6(0)^..^^ . . . ar>* (4.2) 

then defining ip by 
gives 



which proves Proposition [O . 

Thus we have reduced the proof of Proposition f4.1| (hence that of Theorem 
to solving (|4.2|) . Denote by .fT the subspace of M d (g> S' fc+1 (M d ) of tensors 



bji-jk+2 sucn ^ na * 
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Here, and below, S ,fe+1 (IR o! ) denotes the space of totally symmetric tensors of 
order k + 1. Define a map T: K -* S k {R d ) by 



T[bijji...j k ) b. 



Recall the definition of the complete trace from equation ( J3.ll ). To show 



that (|4.2|) can be solved for tp 1 ^ Jfc+i , we use ( p. 10|) and the following Lemma. 



Lemma 4.1. An element b e S (M. ) m t/je image of T if and only if 
ctr fc (6) = 0. 

Proof. Calculation shows that ctr^ o T = 0. To show that the kernel of ctr^ 
is contained in the image of T, we shall use a decomposition ( |4.3|) of 5" (R ). 
For any symmetric product S 9 (M. d ) we define S^R^) as the set of trace-free 
tensors, that is, 

S^(R d ) = {Z n .., q e S%R d )\8^Z m ___ jq _ 2 = 0}. 

If k — q — 2m then we can find a copy of S*^(R d ) in S k (M d ) as the image of 
the map 

W h ...i q i ► 5 Ulj2 . . . S hm _ lhm Wj 2m+1 ... jk ). 

A calculation (or the corresponding decomposition of 0(<i)-modules given in 
§19.5 and Lemma 17.15 of Fulton & Harris (1991)) shows that S k (R d ) is the 
direct sum of all such images. Thus we have 

S k (R d ) = S [k] (R d ) © S [k - 2] (R d ) © . . . , (4.3) 

ending in S^R^) or 5 , ^(R a! ), depending on whether k is odd or even. If 
b G S k (M> d ) then the complete trace of b is either zero (if k is odd) or a 
non-zero multiple of the component of b in ^'(R^) (if k is even). 

It suffices to show that every tensor in the irreducible component (R d ) 
for < q < k is in the image of the map T. Let 

be such a tensor. Then consider 

Xii 1 ...i k + 1 — 5i(^i 1 . . . 0~i 2m i 2m + 1 Zi 2m + 2 ...i k + 1 ) — 5(^42 • • • ^2m + l«2m + 2 ^2m + 3---«fc + l)i- 

This clearly satisfies 

A calculation shows that 

,_ d+q-2 
J- \^n 1 ...i k+1 ) — ^ ^ ^ii...ifc- 

Since d > 1, this proves the result. □ 
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5 The one-dimensional case 



^(x) = i))x j + -M k x j x k + 



We have already noted that Theorem cannot be true in the case that 
d = 1. To be able to determine when there is a ip such that ip(f,b) = (g,c) 
we must either restrict the possible (/, b) or define extra invariants. We 
consider both possibilities in this section. 

The phylon group V(d) has two connected components. If ip is an element 
of V(d) with components 

1 
2! 

then the different connected components of V(d) are determined by the sign 
of the determinant of the matrix (ipj). Denote by Vo(d) the connected 
component of V(d) containing the identity. This is the component where 
det(V>j) > 0. 

In the case that d = 1 we shall write a general element ip G V(l) using 
different notation as 

1p(x) = ijjiX + -^1p2X 2 + 

Then ip is in Vq(1) if ipi > 0. The important fact about this group for our 
purposes is 

Theorem 5.1. If f and g are elements of Ai(l) then there is a unique ip 6 
Vo(l) such that ipf = g. 

Proof. Morse's Lemma shows that ip exists. To show that ip is unique it is 
enough to show that if <pq = q for some <fi G Vq{1) then 0=1. If <pq = q then 
4> preserves length and, as we have noted previously in one dimension, this 
means that <fi(x) = ±x for all x. For to be in the connected component of 
the identity we must have 0(x) = x, so that = 1. □ 

Notice that in the one-dimensional case it follows from ( |3.10| ) and (|3.11|) 
that for any b G V(l) the invariant Aj(6) is a non-zero multiple of the ith 
coefficient in the Taylor expansion 

b(x) = b + bix + 7y b 2 x 2 + . . . , 

if i is even and is zero if i is odd. 

The definitions of even and odd functions extend also to power series. 
We define a power series / to be even if f{—x) = f(x) for all x and odd if 
f{—x) = —f(x) for all x. If g is any power series we define 

1 1 

9 + (x) = -(g(x) + g(-x)) and g (x) = -(g(x) - g(-x)). 
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Then g(x) = g + {x) + g~(x) is the unique expansion of g(x) as a sum of 
an even power series and an odd power series. An even (odd) power series 
has only even (odd) terms in its Taylor expansion. Consider the subgroup 
-Fodd(l) °f -Fb(l) consisting of all ip which are odd. Then we have 

Theorem 5.2. Suppose that (f,b) and (g,c) are elements of M.(l) x V(l). 
with f and g even and that Aj(/, b) = Ai(g, c) for all i = 0, 1, 2, ... . Then 
there is a if> G P <m(1) such that ipf = g and ipb + = c + . 

Proof. From Theorem |5.1| we have that there is a unique ipi G Vq(1) sat- 
isfying ipif = q. Define ip\(x) = —ipi(-x). Then ipif(x) = f(ipi l {. x )) — 
f(— r tpi 1 (—x)) = f(ipi 1 (—x)) = q(—x) = q(x). Hence by uniqueness ?pi = ipi, 
so that ip\ G Podd(l)- Similarly there is a ip2 G P dd(l) such that ip2g = q. 

By invar iance we deduce that Aj^ife) = Ai(i/j 2 c) and hence that ipib and 
ip 2 c have the same even Taylor coefficients. Thus {ip\b) + = (ip2 c ) + - Consider 
now a b of the form b(x) = x k . Letting dip\ denote the derivative of tp\ we 
have 

so that 

(^b)(-x) = {-Ifi^b)^)- 
It follows that for any b, because b + is even, we must have that ip\b + is even 
and similarly ip\b~ is odd, so that 

ipib = ipib + + ipib~ 

is the decomposition of ip\b into even and odd functions and, in particular, 
(i[>ib) + = if)\b + . Similarly {ip2c) + = ip 2 c + so we have ipib + = ip 2 c + - Letting 
ip = ?/^~Vi gives the desired result. □ 

Consider now the case of general (/, b) and (g, c) with A,(/, b) = Aj(g, c) 
for all i — 0, 1, 2, . . . We proceed as in the proof of Theorem |5.2j . We can find 
ipi and ip2 such that ipif — q and tp 2 g = q. We then have Aj(-0i6) = Aj(-02c) 
for all i = 0,1,2,... and hence the even Taylor coefficients of ip\b and ip 2 c 
equal but we have no information about the odd coefficients. If (/, b) is in 
A4(l) x V(l) we define Xi(f,b) to be the zth Taylor coefficient of xpb, where 
i> G 'Podd(l) an d ipf = q. Then we have 

Theorem 5.3. If(f,b) and (g,c) are elements of A4(l) x V(l) then 

Xi(f, b) = Xi(g, c) for all i = 0, 1, 2, . . . 
if and only if there is a ip in Vo(l) such that 

(f,b)=i>(g,c). 
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Proof. The proof follows from the discussion we have already given that 
motivated the definition of the Aj. □ 



We have no explicit formula for the calculation of the \i(f, b). Let us 
show that an algorithm is possible. This is clear, except for finding the ip 
that solves ipf = q. This equation is equivalent to f(ip~ 1 (x)) = x 2 and hence 
f(x) = ip(x) 2 . If we let 

f(x) = ^f 2 x 2 + ^f 3 x 3 + ... 



and 



fax) = fax + Ty^x 2 + ■ ■ ■ 



then the equation f(x) = if)(x) 2 is equivalent to an infinite system of equa- 
tions 

\h = ft 

^7/3 = fafa , 
3! 5.1 
111 



which we can solve recursively for fa , fa > ■ ■ ■ ■ 



Appendix A: Moment generating functions 

We include here some standard results on moment generating functions which 
are needed to derive formula ( |2.1| ). The normal distribution with mean and 
variance matrix E = (E y ) on M. d has probability density function 

<K*) = (2 7 r)^(detS)V2 ex P^^^ V >' ^ 

where we assume the summation convention and the matrix is the 

inverse of the matrix (E 1 - 7 ). If / is a function, we define its expectation with 
respect to this distribution by 

E[f(x)} = [ f{x)(j){x)dx 1 . . . dx d . 
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The moment generating function M(t), where t = (ti, . . . , to), is defined 

by 

M(t) = ElexpiUx*)}. (A.2) 

By completing the square in ( |A.2| ) and using the fact that the integral of 
( A.l ) over all of of R is 1, we can show that 

M(t) =exp(^H t t j ). (A.3) 

Repeated differentiation of (|A.3|) followed by evaluation at t — shows 
that 

Eix^X^ X ik ) = ^ ^ ^fclfc2^fc3fc4 Y^r-lkr 

where the sum is over all partitions of i\, . . . , i r into subsets of size 2. Hence 
if r is odd this expectation is zero. It is straightforward now to calculate the 
formula in (12.11). 
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